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Abstract 
>■ 

o 

0^ ' We obtain general formulae for the plus- and cross- polarized waveforms of 

f—^ ■ gravitational radiation emitted by a cosmic string loop in transverse, traceless 

(synchronous, harmonic) gauge. These equations are then specialized to the 

case of piecewise linear loops, and it is shown that the general waveform for 

1 \ such a loop is a piecewise linear function. We give several simple examples of 

wo. 



the waveforms from such loops. We also discuss the relation between the grav- 
itational radiation by a smooth loop and by a piecewise linear approximation 
to it. 
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I. INTRODUCTION 

Cosmic strings are one dimensional topological defects that may have formed if the 
vacuum underwent a phase transition at very early times breaking a local U(l) symmetry 
|H|-Q. The resulting network of strings is of cosmological interest if the strings have a large 
enough mass per unit length, \i. If G/z/c 2 ~ 10~ 6 , where G is Newton's constant and c is 
the speed of light (i.e. \x ~ 10 22 g/cm) then cosmic strings may be massive enough to have 
provided the density perturbations necessary to produce the large scale structure we observe 
in the Universe today and could explain the pattern of anisotropics observed in the cosmic 
microwave background ||. 

The main constraints on \x come from observational bounds on the amount of gravita- 
tional background radiation emitted by cosmic string loops ( @J§,0] and references therein). 
A loop of cosmic string is formed when two sections of a long string (a string with length 
greater than the horizon length) meet and intercommute. Once formed, loops begin to oscil- 
late under their own tension, undergoing a process of self-intersection (fragmentation) and 
eventually creating a family of non-self-intersecting oscillating loops. The gravitational radi- 
ation emitted by each loop as it oscillates contributes to the total background gravitational 
radiation. 

In a pair of papers, we introduced and tested a new method for calculating the rates at 
which energy and momentum are radiated by cosmic strings @||. Our investigation found 
that many of the published radiation rates were numerically inaccurate (typically too low 
by a factor of two). Remarkably, we also found a lower bound (in the center-of-mass frame) 
for the rate of gravitational radiation from a cosmic string loop [K| . Our method involved 
the use of piecewise linear cosmic strings. In this paper we wish to provide greater insight 
into the behaviour of such loops and, in particular, how they approximate smooth loops by 
examining the waveforms of the gravitational waveforms of such loops. 

It has long been known |||lT) that the first generation of ground-based interferomet- 
ric gravitational-wave detectors (for example, LIGO-I) will not be able to detect the 
gravitational-wave stochastic background produced by a network of cosmic strings in the 
universe. The amplitude of this background is too weak to be detectable, except by a future 
generation of more advanced instruments. However, a recent paper by Damour and Vilenkin 
[|T^] has shown that the non-Gaussian bursts of radiation produced by cusps on the closest 
loops of strings would be a detectable LIGO-I source. While the specific examples studied 
here do not include these types of cusps the general method developed can be applied to 
such loops. 



Our space-time conventions follow those of Misner, Thorne and Wheeler [|13| so that 
i]^ v = diag(— 1, 1, 1, l) M y. We also set h — c — 1, but we leave G explicit. 

II. GENERAL THEORY 

In the center-of-mass frame, a cosmic string loop is specified by the 3-vector position 
x(£, a) of the string as a function of two variables: time t and a space-like parameter a 
that runs from to L. (The total energy of the loop is fiL.) When the gravitational back- 
reaction is neglected, (a good approximation if Gjj 2 <K 1), the string loop satisfies equations 
of motion whose most general solution in the center-of-mass frame is 



x(t,or) = i[a(ti)+b(t;)]. (2.1) 

where 

u = t + a, v = t- a. (2.2) 

Here a(-u) = a(w + L) and b(i>) = b(v + L) are a pair of periodic functions, satisfying the 
"gauge condition" |a'(u)| = |b'(v)| = 1, where ' denotes differentiation with respect to the 
function's argument. Because the functions a and b are periodic in their arguments, the 
string loop is periodic in time. The period of the loop is L/2 since 

x(t + -, a + -) = -[a(t + a + L) + b(t - a)] = -[a(t + a) + b(t - a)] = x(t, a). (2.3) 

With our choice of coordinates and gauge, the energy-momentum tensor T^ v for the 
string loop is given by 

T""(t,y)=// ff dudvG^(u,v)5(t-±(u + v))5 {3) (y-x(u,v)), (2.4) 

0<u-v<2L 

where G^ v is defined by 

G^(u, v) = d u x»d v x u + d v x»d u x\ (2.5) 

with x^ = (t, x(t, a)). In terms of a and b, 

G 00 = i G 0i = \[a[ + b% G» = \[a% + a'b>], (2.6) 

and the trace is 

G A A = |[-l+a / -b']. (2.7) 

Alternatively we may introduce the four-vectors A^{u) = (u, a(-u)) and B^(v) = (v, b(v)) so 
that 

qhv = \^b w + B'^A'"). (2.8) 

The "gauge conditions" are satisfied if and only if A ,fl (u) and B /fl (v) are null vectors. 

As a consequence of the time periodicity of the loop the stress tensor can be expressed 
as a Fourier series 

CO 

T ltv (t,y)= £ e w T^K,y), (2.9) 

n=— oo 

where u n = Ann/ L and 

L/2 

2VK,y) = | J dte-^X u {t,y) 

o 

L L 
= -j- I du f dve-^^ u+v) G^(u,v)5 {3) (y -^(u,v)). (2.10) 





The retarded solution for the linear metric perturbation due to this source in harmonic 
gauge is 
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V(*,x)=4G J2 



gy 

l x ~y| 



T^ v (uj n ,y) - IrifvT \(w n ,y) 



p iw„(t-|x-y|) 



(2.H) 



Far from the string loop center-of-mass the dominant behavior is that of an outgoing spher- 
ical wave given by 



V(^x)=4G J2 



oo iu n (t-r) 



d y T^(u n ,y) - ^q^T \{u n ,y) 



Ju) n fl-y 



(2.12) 



where r = \x\ and 17 = x/r is a unit vector pointing away from the source. Inserting Eq. 
( p.lOj) into Eq. ( |2.12 ) we find the field far from a cosmic string loop is 



, , , 8Gu ~ e i^(t~r) j j 
/i MI/ (£, xj = — — 2_^ I du I dv 



L n=-oo r 



G^(u,v) - pi/jvG x (u,v 







(2.13) 



The n = term in this sum corresponds to the static field 

L L 



/ l st f ic (t,x) = ^ / (I// / dr 



rL 



G^v(u,v) - \r]^ u G X x{u,v) 







(2.14) 



2G/iL 



{jiya, + 2t fM i u ) 



2GM 



IMU1 



(2-15) 



as appropriate to a object with mass M as may be seen by comparison with the Schwarzschild 



metric in isotropic coordinates (see, for example, Eq. (31.22) of Ref. [0)- We denote the 
radiative part of the field by 



trad -L _ 7 static 



We may rewrite Eq. fl2.13|) as 

oo 



(2.16) 



(2.17) 



where k^ = (1, ft) is a null vector in the direction of propagation and 

L L 



» 



fiu 



8GM 
rL 2 



du / dv 



G^(u,v) - \rj^G x x {u,v) 



s iicj n [fe M ^(u)+fe M BM(«)] 



(2.18) 







are polarization tensors. From Eq. (|2.8| ), it is clear that the polarization tensors may be 
written in terms of the fundamental integrals 



1 L r 
jWm = - du A'»{u) e^ k » A ^ u \ (2.19) 

o 



and 

L 







In terms of these integrals 



2GM 
eoo = [Wo + I J] (2.21a) 



2GM r , 

e 0i = [I Ji + J Q Ii] (2.21b) 



c 



_ 2GM 
r 



{[IiJj + Jilj] + Sij [IqJq - I-J]} , (2.21c) 



where we have dropped the superscript (n) for clarity. 

The harmonic gauge condition requires that the polarization tensors satisfy k^e^ = 
^k^e^. This is easily verified by noting that 1° = fM and J° = fl-J. These equations 
follow from the identity 

L L 

f du hj^iu) e-^ nKAV{u) = — I ' du^- e-^nk.A^u) = Q ^ ? 2 ^ 

J to J du 



which is a consequence of periodicity, and the corresponding equation for B^. The harmonic 
gauge condition does not determine the gauge completely and we are left with the freedom 
to make transformations of the form 

G fxv — 6;ui/ T Kpjiy + K V E ^. yZ.Zoj 

If we make the choice 

£o = — [Wo + I-J] (2.24a) 

r 

and 

e t = — {[/ Jo + I-J] a + 2 [I J t + J /,]} (2.24b) 

r 

then 

eV = 0. (2.25) 

The spatial components are given by 



2GM 
e'ij = {[hJj + IjJi] + 5 lj [hJ - I J] + ^QjlIoJo + I- J] 

+ IolJitlj + Q l J j ] + Jollify + SUIj]}, (2.26) 

these satisfy the gauge conditions 

e'/ = e'> = 0, (2.27) 

and 

e'/Qi = 0. (2.28) 

If we perform a spatial rotation to coordinates (a/, j/', z') where Cl points along the z'-axis 
then we can write 

e+ e x 0' 
e' iY = I e x -e+ I , (2.29) 



where 



and 







2GM 
e+ = [IyJy - I V J V ] (2.30a) 



2GM 

e x = [I V J 2 , + h'Jvl (2.30b) 

r 



define two modes of linear polarization. 

In terms of the original basis we can write 



e+ = (cos2^v4+ + sin2^v4 x ) (2.31a) 



and 



2GM 

e x = (-sin2^A + + cos2^A x ) (2.31b) 

r 

with 

A+ = I (cos 2 cos 2 — sin 2 (j))I\Ji + (cos 2 9 sin 2 — cos 2 0) J2 J2 + sin 2 6*73 J3 
+ (1 + cos 2 0) cos</>sin0(/ 2 </i + /1J2) — sin#cos0cos</>(i3Ji + I\J 3 ) 
- sin cos 9 sin 0(/ 2 J 3 + J 3 J 2 ) } (2.32a) 

A x = J —2 cos 6* sin cos §I\ J\ + 2cos6*sin</>cos0/ 2 ^2 + cos #(2 cos 2 — l)(J 2 Ji + I1J2) 

+ sin9sm<j){I 3 J 1 + I 1 J 3 ) - sin 9 cos <p{I 2 J 3 + I3J2)} (2.32b) 

where 9, <ft and ip are the Euler angles defining the orientation of the frame (x', y', z') relative 



to the original frame (our conventions follow those of Ref. [15]). The corresponding linearly 
polarized waveforms are then defined by 



X 



h™f x {t,x) = ^{e i47rn(t_r)/L e^ x + e-^ n(t - r)/L e^;}. (2.33) 

n=\ 

Recall that /i rad is obtained from the full metric perturbation h by dropping the n = term 
which corresponds to the static (non-radiative) part of the field. 
The power emitted to infinity per solid angle may be written as 

HP r 2 

__ = ii m -(h a0t h a0 r ) (2.34) 

rij/2 oo 

= ^E^{(l4 n) | 2 ) + (l4 n) | 2 )}. (2.35) 

27r n=l 



III. EXAMPLES 

For convenience we shall now set the length of the loop L = 1, and take ■0 = 0. 

A. Piecewise linear loops 

These are loops for which the functions a(-u) and b(t>) are piecewise linear functions. 
The functions a(w) and b(t> ) may be pictured as a pair of closed loops which consist of 
joined straight segments. The segments join together at kinks where a'(-u) and b'(v) are 
discontinuous. 

Following the notation of Ref. || we take the a- and 6-loops to have N a and iV b linear 
segments, respectively. The coordinate u on the a-loop is chosen to take the value zero at 
one of the kinks and increases along the loop. The kinks are labeled by the index i where 
i = 0,1, . . . ,N a — 1. The value of u at the zth kink is denoted by U{ and without loss of 
generality we set Uq = 0. The segments on the loop are also labeled by i, with the ith 
segment being the one lying between the ith and (i + l)th kink. The kink at u = u^ a is the 
same as the first kink at u = uq = but, even though uq and UN a are at the same position 
on the loop, uq = while UN a = 1- The loop is extended to all values of u by periodicity 
(with period 1). We denote a, = a(wj), and the constant unit vector tangent to the ith 
segment by a^. Then we have 

a(w) = a, + a-(w - -Uj) for u E [ui, Ui+i], (3.1) 

and for consistency 

a i+ i = aj + &' i (u i+1 -u^. (3.2) 

We have corresponding definitions for the b-loop and we follow the convention of Ref. by 
labeling the kinks by the index j. 

It is now elementary to calculate that, for n^O, 

1 Ara_1 a' 

j(") _ V^ i ( -2irin(u i+1 -Sl-a i+1 ) _ -2mn{ui— O-Bi) \ 

" 2mn ^ 1 - £1.^ I J 

i N ^ { a' a' 



V^ J gj i-l I e -2irin(tti-n-ai) /g g\ 



2-im ^ [1-fl^ l-n-aj.! 



with a similar equation for J. If we insert these expressions into Eq. ( |2.32| ) and then into 
Eq. (|2.33|) the sum over n for h x j& L consists of terms of the form 

OO 1 

53 — cos 2vrn[2(t - r) + (m - fi-a,) + ( Vj - fl-bj)] (3.4) 

n=l n 

which may be performed exactly using the identity 

OO -| 

y^ — cos 2irnx = tt 2 (x 2 — X + g) l£ [0,1). (3.5) 

n=l n 

This function is extended to other values by periodicity, for example, for a; G [1,2) we 
merely replace x by i - 1 in Eq. ( |3.5| ). Such transformations leave the coefficient of x 2 
unchanged and can only change the coefficient of x by a multiple of 2. As a result when 
the sum in Eq. ( |3.3| ) is performed for the coefficient of x 2 the sum telescopes and gives 
zero. Thus, the waveform of a piecewise linear loop will be a piecewise linear function. In 
addition, considering the coefficient of x all slopes of the waveform must be a multiple of 
some fundamental slope. The slope only changes when a (4-dimensional) kink crosses the 
past light cone of the observer at (t, x). These properties are illustrated in the examples 
below. 



B. Garfinkle-Vachaspati loops 



As our first set of loops we study the loops considered by Garfinkle and Vachaspati JIB . 
The vectors a(w) and b(v) lie in a plane and make a constant angle 2a with each other 
where a G (0, tt/2). To be specific, we may take a(-u) and b(y) to be given by 

w(cosai + sinaj), u G (0, |) 

a(-u) = { (3.6a) 

(1 — m) (cos ai + sinaj), u G (|, 1) 



b(u) = <^ 



f (cos ai — sinaj), u G (0, |) 

k (1 — v) (cos ai — sinaj), wg(|,1) 
It is then straightforward to calculate that, for n^O, 



(3.6b) 



piirnyl— sin6>cos(</)— «)) i 

I (n) = - — r, r^-7 ^7-7 r- (cos ai + sinaj) (3.7a) 

«7m(l — sin o cos" : (0 — a)) 

p iirn\l—sin9 cos(<j>+a)j i 

J (n) = - — r-j- j-— — -(cosai - sinaj) (3.7b) 

«7m(l — sm ^cos z (0 + a)J 



and correspondingly 



A 



(«) 



sin 2 9 cos(0 + a) cos(0 — a) — cos 1 



x 



1 — sin 2 9cos 2 ((f) — a))(l — sin 2 9 cos 2 ((f) + a)) 

' inn[l—sin8cos(<f)—a)j i \ f innyl— sin0cos(</>+a) ) _ 



n 2 TT 2 



A 



(n) 



l)(e« 

cos 6 1 sin 20 



x 



1 — sin 2 #cos 2 (0 — a))(l — sin 2 #cos 2 (0 + a)) 

-i \ / iirny i— sin t/ cos (^H-an _ i\ 



' J7rn(l— sin6cos(</>— a)) -i \ / innyl— sin0cos(</)+a) ) 



n 2 7T 2 



(3.8a) 



(3.8b) 



As described above, the sum over n in Eq. ( |2.33|) may be performed explicitly to yield a 
piecewise linear function. For example, G [0, vr/2), h + is given explicitly by 

sin 2 9 cos(0 + a) cos(d> — a) — cos 26 



h, 



x 



2GM 

r (1 - sin 2 6 l cos 2 (0- a))(l - sin 2 #cos 2 (0 + a)) 
' (1 — sin6*cos(0 — a))(l — sin6'cos(0 + a)), 

< (t — r) < | sin 6'(cos(0 — a) + cos(0 + a)) 

-8(t-r) + (l + sin6>cos(0 + a))(l + sin6>cos(0-a)), 

^ sin 9 (cos (cf) — a) + cos(0 + a)) < (t — r) < |(1 + sin 6* cos(cf) + a)) 



A 



sin^ cos(0 — a))(l + sin^ cos(0 + a)), 

i(l + sin6 l cos(0 + a)) < (t-r) < \(l + sin6 l cos(0 - a)) 



Kt 



2) + (1 — sin6'cos( 



— a))(l — sin6'cos( 
7(1 + sin6'cos( 



+ «)), 
- «)) < (* 



r)<|. 



(3.9) 



and the waveforms are periodic in t with period |. The intervals are ordered in the given 
way for our choice of G [0, vr/2). h x is obtained simply by replacing the prefactor by that 
appropriate to A x as is clear from Eq. ( |3.8| ). To obtain the waveforms for other angles we 
may note that the transformation <p — > <fi + tc is equivalent to changing the sign of (t — r), 
while the transformation — > 7r — is equivalent to changing the sign of (t — r) and sign in 
front of the sin 20 term in the prefactor in h x . 

Note that the apparent singularity in the waveforms in the plane of the loop (9 = 7r/2) 
at = ±« and = it ± a is spurious. This may be seen by noting that the waveform is 
bounded by the two constant sections of the piecewise linear curve which take on a value 
which tends to zero in this limit. In fact, the numerator of the prefactor also vanishes in this 
limit which ensures that the amplitude tends to zero at these points and hence that even 
the time derivatives (which determine the power) are finite. Along the axis 9 = 0, Eq. ( |3.9| ) 
reduces to 



h +/x (9 = 0) 



2GM 



cos 20/ sin 20 



l(t-r)+l, 0<t-r< 



\(t-r-\) + l 



(3.10) 



4 — l ' ^ 2" 



Waveforms for various angles a plotted in Fig. [I] for the case of a = it/A, corresponding 
to two lines at right angles. This is the configuration which radiates minimum gravitational 
radiation for this class of loops, P = 64 In 2 Gjj 2 w 44-3614 G[i 2 . 
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C. Plane-line loops 



As our next set of examples we study the set of loops in which a(w) lies along the z-axis 
and b(v) is always in the x-y plane. This class of loops was studied by us in Ref. fl?| 



where we gave an analytic result for the power lost in gravitational radiation by such loops. 
Explicitly a(w) is given by 



uk 



a(u) 



«e(o ,|) 



;i-«)k ue (|,1) 



It follows that 



A-KTlll — COS 0) 1 

tW - - _k 

■ 2 n 



Also A n) 



0, so we have 



(3.11) 



(3.12) 



A, 



sin 6 cos #(cos <fi J\ + sin 4>J2)h 



A x = sin^(sin0Ji — coscfi^h- 



(3.13a) 
(3.13b) 



It follows immediately that the waveforms vanish along the z-axis. 

In Ref. [|17| we proved that the minimum gravitational radiation emitted by any loop in 
this class is given by taking the b-loop to be a circle: 



b(v) = — (cos(27ru)i + sin(27Tf )j). 

27T 

The power emitted in gravitational radiation by this loop is 

27T, _ - 

P = 16 / [ C ° SX) dx G/i 2 « 39-0025 G/x 2 . 
j a; 

o 



J*™) may be determined explicitly as 



A 



(n) 



J, 



(n) 



3 <(«+i)(«-|)j n+1 ( n8 i n 0) + e *(™-i)(*-f)j n _ 1 ( ns i n 0) 



This gives the equivalent forms 



.4 



(n) 



sin(7msin 2 (#/2)) COS6 1 



(3.14) 



(3.15) 



(3.16a) 
(3.16b) 



-2 



sin# 
sin(7msin 2 (6>/2)) cos(9) 
sin 2 # 



[J n+ i(nsin0) + J n _i(nsim 
1 



3 i7rnsin 2 (e/2)+m(</)-7r/2) 



7rn 



J ra (nsinfcM — < 
Tin 



iruf>—in(ir/2) cos ( 



(3.17) 



and 
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4 n) = / n(7mSin ^ /2)) [J n+1 (nsm9) - J^nsine)} ± e **»*?V/*W*-*/*) 
sin 9 im 

= 2 ^nsmy/2)) g) J_ e ^_ w(7r/2)OTs9 g 

sinfl ' im ' 

The corresponding waveforms for various choices of 9 are plotted in Figs. and || (As the 
system simply rotates cylindrically with time the choice of <fr is irrelevant, corresponding 
simply to a shift in t — r.) 

In the plane of the b-loop h + vanishes so that the wave becomes linearly polarized. 
On the other hand, as we approach the axis 9 = the fundamental mode (n = 1 term) 
dominates and we have 

h+ ~ sin#sin(47r(£-r) + 0) (3.19) 



and 



h x sin#cos(47r(£-r)+0). (3.20) 



Thus the wave approaches circular polarization but its amplitude vanishes as sin 9. 

As in Ref. [O] we may also consider the case where the b-loop forms a regular iV-sided 
polygon. In Figs. f| and || we compare the waveform for the circle with that for a regular 
hexagon for which P = (40 In 5 — 32 In 2) Gfj, 2 ~ 42-1968 Gfi 2 . As mentioned above a change 
in for the circle-line loop corresponds simply to a shift in t, however, this is no longer the 
case for the polygon for which the waveform will only repeat every 2n/N. Hence in Figs. |] 
and |5] we include hexagon-line waveform for both = and (f) = 5ir/6 (this choice was made 
simply to disentangle the two graphs as far as possible). It is remarkable that even for such 
a crude approximation to the circle as a hexagon, the waveform of the hexagon-line loop 
provides remarkably good piecewise linear approximations to the circle-line waveforms. 

IV. CONCLUSION 

Given the remarkable agreement of the waveforms it is of interest to compare the 'in- 
stantaneous power' defined by 

_ GM 2 ~ 2 (w) 2 

Z7r n=l 

in the different polarizations. While this quantity is not gauge invariant its time average 
is and gives the total power radiated in each polarization. By comparing the function for 
the polygon-line loops with the circle-line loop we can certainly see that their time averages 
agree well. As the waveform for a piecewise linear loop is a piecewise linear function, the 
instantaneous power, which is the square of its derivative, will be piecewise constant. For 
example, in Fig. || we compare the 'instantaneous power' in the plus-polarization between 
the circle-line loop and a regular 24-sided polygon-line loop. The very close agreement 
between the two curves provides further evidence for the validity of the piecewise linear 
approximation of string loops used by |§ . 
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FIGURES 
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FIG. 1. Plus-polarized waveforms for the Garfinkle-Vachaspati with a = ir/4. The solid line 
corresponds to the wave travelling up the z-axis (6 = 0). The dotted line corresponds to a direction 
at elevation 9 = tt/3 along (f> = 0. The dashed line corresponds to a wave travelling in the plane 
of the loop 9 = 7r/2) at angle </> = tt/5. The cross-polarized waveforms differ only in that there 
amplitude has a different dependece on 9 and (/>. 
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FIG. 2. Plus-polarized waveforms for the circle-line loop. Plotted is h + /svn9 for 9 = tt/2 
(solid), 9 = ir/4 (dashed) and 9 = tt/20 (dotted). The choice of scaling is chosen on the basis of 
the asymptotic form Eq. ( |3.19| ). 
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FIG. 3. Cross-polarized waveforms for the circle-line loop. Plotted is /i x /sin# for 9 = tt/2 
(solid), 9 = 7r/4 (dashed) and 9 = 7r/20 (dotted). The choice of scaling is chosen on the basis of 
the asymptotic form Eq. (|3.20| ). 
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FIG. 4. The solid lines are the plus-polarized waveforms for the hexagon-line loop with 9 = 7r/4 
and with (f> = (left through) and <p = 5ir/6 (right through). The dotted lines are the corresponding 
waveforms for the circle-line loop. 




FIG. 5. The solid lines are the cross-polarized waveforms for the hexagon-line loop with 8 = 7r/4 
and with <fi = (left peak) and 4> = 5ir/6 (right peak). The dotted lines are the corresponding 
waveforms for the circle-line loop. 
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FIG. 6. Comparison of the 'instantaneous power' in plus-polarized waves for the 24-sided poly- 
gon-line loop (dotted line) and for the circle-line loop (solid line) with 9 = 7r/4 and = 0. 
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